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We introduce the Wigner functional representing a quantum field in terms of the field amplitudes 
and their conjugate momenta. The equation of motion for the functional of a scalar field point out 
the relevance of solutions of the classical field equations to the time evolution of the quantum field. 
We discuss the field in thermodynamical equilibrium and find the explicit solution of the equations 
of motion for the so-called "roll-over" phase transition. Finally, we briefly discuss the approximate 
methods for the evaluation of the Wigner functional that may be used to numerically simulate the 
initial value problem. 



I. INTRODUCTION 

Dynamical processes in relativistic quantum fields are usually described in terms of the particle excitations of 



CO 

> 

' these fields. However, there are situations where more appropriate degrees of freedom are rather the field amplitudes 
themselves and their conjugate momenta. This happens when as a first approximation the dynamics is well described 
• by the classical field equations. A well known example here is the evolution of Higgs fields in the early universe [0,0, 
which acquire a finite vacuum expectation value when the temperature falls below the critical temperature for the 
symmetry breaking phase transition. The mechanism which drives the Higgs field from the symmetric phase, where 
the vacuum expectation value vanishes, to the asymmetric phase with a non-zero expectation value is essentially 
classical. Such phase transitions have attracted much attention in the context of inflationary cosmological models 
H, models of galaxy formation Q, the creation of the cosmic baryon number asymmetry during the electroweak 
O i| phase transition ||, and most recently due to the possibility of creating the so-called misaligned chiral condensates 
in ultrarelativistic heavy-ion collisions ||. 

All these situations have in common that one must study the evolution of a quantum field far off thermal equilibrium, 
when the initial state of the quantum field is specified. Most theoretical approaches to this initial value problem for 
quantum fields have been based on dynamical equations for the field expectation values and the Gaussian fluctuations 
around those assuming that fluctuations around the most probable path remain small. While this assumption 

■ may be correct in certain instances, it is a well known fact that fluctuations usually do not remain small in dynamical 
phase transitions, where domain formation and clustering are regularly occurring phenomena. 

In this paper we present a formalism which appears to be very convenient to study the evolution of quantum fields 
from an initial state and which goes beyond the classical description. The central object of our approach is the 
Wigner function, Wigner functional in fact, which is the quantum analog of the classical distribution function in a 
phase space. Here the phase space, however, is spanned not by the particle coordinates and momenta, but by the field 
amplitudes and their conjugate momenta. The great advantage of the Wigner function formalism is that, while being 
fully quantum mechanical, it remains close in its spirit to the classical description, and that the classical limit can 
always be easily obtained. On the other hand, the simultaneous presence of many "quasiclassical" field configurations 
can be described with ease, and the fluctuations around these do not have to remain Gaussian. 

An additional motivation for our work derives from the fact that there are many interesting investigations of 
nonlinear classical field equations in Minkowski space. One can mention here studies of the chaotic behaviour of 
classical, in particular, Yang-Mills fields ||, searches for exact solutions of the equations of motion (see e.g. p0| for 
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recent developments), or very specific predictions concerning multiparticle production in high-energy collisions 
Whereas the relevance of classical solutions in euclidean space is safely grounded because they minimize the euclidean 
action, as in the case of instantons |J, no such a argument exists for solutions in Minkowski space. Attempts to 
understand the quantum nature of classical fields have been based on the coherent state representation of Fock space 
lyl or semiclassical methods (juj). For recent work in this direction, see JTif . 

In order to avoid unnecessary complications we restrict considerations in this paper to scalar quantum fields with 
quartic self-interaction. After discussing the definition of the Wigner functional we show that the equations of motion 
have the familiar form of a transport equation in phase space with quantum corrections. We then derive the Wigner 
functional for a free scalar field in thermal equilibrium. We discuss how the classical phase space distribution is 
obtained in the high temperature limit, and we obtain the two-point correlation function. We analyze the "rollover" 
of the scalar field in a second-order phase transition, where the power of the Wigner functional approach becomes 
apparent, allowing for the evolution of the quantum field along many simultaneous classical paths. We briefly discuss 
the mean field and semiclassical approximation, and finally we conclude with a suggestion how to obtain numerical 
solutions to the initial value problem for quantum fields in more complex situations. 



II. WIGNER FUNCTIONAL 



We consider a real scalar quantum field with the standard Lagrangian density 



£{t,x) = ^d^{t,x)d„${t,x) - ^m 2 &(t 7 x)+£ I (t,x) . 



(1) 



The interaction Lagrangian £j is assumed to be a polynomial in $ such as £j — — ^-$ 4 . We use hats to distinguish 
operators from corresponding c— numbers. For simplicity most formulas are written for the 1 + 1 dimensional field 
theory, but the generalization to higher dimensions is straightforward. 

Before defining the Wigner functional of the field, let us write down the well known definition of the single particle 
Wigner function as fl^-|l8|] 



W(q,p;t) 



due-^(q+ l -u\p(t) \q- 



1 



(2) 



where p and q are the particle momentum and position, respectively, while p is the time dependent density matrix 
operator in the Schrodinger picture. In analogy to (||) we define the Wigner functional as 



W[$(x),H.(x);t] = I Vcf>'(x) exp 

1 



-i I dxU{x)<p'(x) 
1 



x{9(x) + -4/(x)\p(t) \$(x)--<j/(x)) 



(3) 



It sometimes appears easier to compute the Wigner functional of the fields in momentum space. However, we then 
face a little complication. The field <j>(x) is real while is complex, but the real and imaginary components of 
are not independent from each other because of the constraint $(— p) — &*(p). Thus, we adopt the following 
procedure. The real and imaginary parts of &(p) are treated as independent variables, but p £ (0,oo) instead of 
p G (—00,00). The Wigner functional is then defined as 



wm P ),iL(p);t] 



i dp(ir(py(p) + n(p)$'*(p)) 



T>4>' (p) exp 



(4) 



with the functional integrations running over real and imaginary components of 4>' (p) . The transformation from (||) 
to (Ml) involves a nontrivial Jacobian 



det 



6$(x) 
6<f>(p) 



(5) 



However, the determinant does not depend on the fields because the Fourier transformation is a linear mapping; 
consequently, the Jacobian affects only the normalization of the Wigner functional. 
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To clarify the physical meaning of the Wigner functional, we consider the expectation value of an operator ft). 
The expectation value is defined as 

(0($,n)) = |Tr[p(t)0($,iT)] , (6) 

with 

Z = Trp= / X>$(» — W[*,II;t] . (7) 

As is shown in Appendix A, 

(o($,n)) = (o($,n)) = - / o($,n) w[*,n-,t] , (8) 

provided that the noncommuting operators in O are properly symmetrized. Operators corresponding to an asymmetric 
ordering of the operators $ and IT must be explicitly expressed as sums of symmetrized terms. 

Equation (|J) shows that the Wigner functional has the same role in quantum field theory that the density distribu- 
tion in the phase space spanned by <& and II has in the classical field theory. We emphasize, however, that the Wigner 
functional is not always positive definite and thus cannot be interpreted as a probability density in phase space. The 
analytically tractable cases discussed below, when the Wigner functional is of Gaussian form in the variables $ and 
II, are somewhat exceptional in this respect. Of course, the limited classical interpretability of the Wigner functional 
does not diminish its usefulness as a representation of the quantum mechanical density matrix with a simple classical 
limit. This will become obvious in the next section, where we study the relevance of solutions of the classical field 
equations for the time evolution of the Wigner functional. 



III. EQUATION OF MOTION 



The density matrix operator satisfies the equation of motion 



(9) 



In this section we refrain from setting h = 1, because we want to discuss the classical limit of the evolution equation. 
The Hamiltonian H for the real scalar field is 



H = 



~ J dx(ti?{x) + (V$(x)) 2 + m 2 &(x) - C z {x) 



(10) 



where II = — $ denotes the conjugate momentum operator. 

As discussed in Appendix B, one derives from Eq. (ph the following equation of motion of the Wigner functional 



where 



d_ 

dt 



^ -(»'•«-*•«) jab 



W[$,II;t] = 



)d(x) = -~Ci $0 



ih S 



2 5IL(x) J h 



ih s 
~2 5U(x) 



For £j($) = -jf$ 4 we find 



(11) 



(12) 



K z = ■- ( -&{x) 



8 ft* S 3 
SU(x) + 4 [x) 5U 3 (x) 



(13) 



In particular, the interaction term K,\ always terminates for renormalizablc quantum field theories in (3+1) space-time 
dimensions, since their Lagrangians contain at most quartic interaction terms. 
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One sees that the equation ([□]) has the familiar structure of a transport equation. When the higher derivative term 
in the interaction (n3|) is neglected, eq. (Ill]) can be written in the form of the classical Liouville equation 



■tt + I dx 



( 8H 



5H 



W[$,Il;t] = 



(14) 



where H is the classical Hamiltonian. The neglected term is proportional to K , showing that corrections to the 
classical phase space evolution are of second order in h. We will later (in section 7) briefly discuss how the quantum 
corrections to the time evolution of the "classical" Wigner functional could be obtained. 



IV. FREE FIELDS IN THERMAL EQUILIBRIUM 



Let us consider the scalar field in thermodynamical equilibrium. The density operator is then 

B=~e- f>6 



(15) 



where (3 = T _1 is the inverse temperature and Z — Tre~^ ff is the partition function. Since there is no conserved 
charge carried by the real scalar field, there is no chemical potential in Eq. (15). 

It appears easier first to compute the Wigner functional of the fields in the momentum space (^) and then transform 
it into the coordinate space representation (0). The Hamiltonian 



H^f™^ (ft f (p)ft(p) + (P 2 + m 2 )&(pMp)) 



(16) 



can be identified as the sum of independent harmonic oscillator Hamiltonians each representing the mode of momentum 
p. Using this analogy one finds, as shown in Appendix C, the equilibrium Wigner functional as 



W^[$,n] = Cexp 



-i/3 

r 



(17) 



with E(p) — \J p 2 + m 2 , the thermal weight factor 



A,(p) = ^-th^M 



and the divergent normalization factor 



C = exp 



*mth« 

2tt 2 



(18) 



(19) 



which represents the contributions from zero modes. Since the expressions under the momentum integrals in ( |l7j ) and 
(HJ) are even functions of momentum, the integrals are extended from — oo to oo not from to oo, as in the definition 
(]]). The Wigner functional ( |l7j ) is normalized to unity. 
In the high temperature limit, where C, A^(p) — > 1, we reproduce the classical result i.e. 



Wg'[$,n] = exp 



\p J £ (n*(p)U(p)+E 2 (p)^(pMp) 



(20) 



In this case the Wigner functional depends only on the total energy of the system. 
One immediately finds from (|17|) that 



(*t(p)$(p)) 



1 



2 E (p) th^M ' 



(tf (p)ft(p)> = -^i 



2 th 



PE{p) 



(21) 



where (...) is given by Eq. (0). Therefore, 



($t(p)$(p)> (nt(p)n(p)) 



2 th 



i l 

Wip) ~ 2 



(22) 
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showing that the uncertainty principle is built into the Wigner functional. 

Knowing the Wigner functional (Q) one easily finds the Wigner functional (^). As mentioned previously, the 
Jacobian (pi) only modifies the normalization factor, which is infinite anyway. Thus, 



f3 dx dx 1 Tt(x, x') 



with 



where 



Wp[<&{x),TL(x)) = C exp 



H(x, x') = -A(x- x') [IL(x)Il(x') + V4>(x)V$(x') + m 2 ${x)<f>(x')) , 



(23) 



(24) 



(25) 



A/9(x) = f°° ^ e"^ K p {p) . 

One sees that A | a(x) w for m _1 , 3> /? which corresponds to the classical limit. For m — the integral (|2 
can be computed analytically in one dimension JljJ as 



A,(x) = - Incth^i, 



and approximated in the following way 



2 l \x\ 

0-k ,8 



exp 



H?) 



for |x| < /3 
for |x| > /3 



(26) 



(27) 



One easily checks that the Wigner functional (|23|) represents a time-independent solution of the equation of motion 
©■ 



V. EQUILIBRIUM CORRELATION FUNCTION 

As usually, the correlation function is obtained from the generating functional, which for the free fields in equilibrium 
is defined as 



Z[j(x)} = J £>$^ W[^,U] exp (b J dx $(x)j(x) 

X>$^2 exp -/3 J dx dx' (H(x, x') - S(x - x')$(x)j(x') 

where H(x,x') is given by Eq. ( |24| ) and j(x) denotes the external classical current. Due to Eq. 
function can be expressed as 



(*(*)*(!/)) 



1 1 S 2 Z[j] 



^ Z\j] 5j(y)5j(x) 
As shown in Appendix D, the generating functional is 



j(x)=0 



Z[j(x)} =AAexp 
where Af is the normalization constant and 



dx dx' j(x)Q(x — x')j(x') 



0(*)=/£t" 



-ipx 



2tt K {p) [p 2 +m 2 ] ' 
with Kp(p) given in Eq. (|Ig|). Thus, the correlation function ( p9| ) equals 



(28) 

the correlation 
(29) 

(30) 



(31) 
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($(x)<%)> = Q(x - y) . (32) 

The integral ([u]) can be evaluated analytically [[f9| in three dimensions for m = 

1 sh(27r|x|//3) 
^ Xj_ 4tt|x| ch(27r|x|//3) - 1 ' 



4^ for /3«|x|, 
f0r /3>> l x l- 



S«H 4 T ^ ;.JJ (34) 



For a finite mass the integral (|31j) can be calculated in an approximate way and for three dimensions one finds 

0(x) w j^t e- m|x| (35) 
47r|x| 

for to" 1 , |x| 3> /3, which corresponds to the classical limit, and 

^ " I^xj Xl(m ' X|) " e_m|X ' ' (36) 

where the first approximate equality holds for mT x <C f3, while the second one requires additionally mT x <C |x|. 

VI. "ROLLOVER" PHASE TRANSITION 

Let us consider the scalar field in variable heat bath, undergoing a phase transition with falling temperature. 
The mass squared of the field, effectively being temperature dependent, is initially assumed to be positive. We are 
interested in the evolution of the quantum field with spatially homogeneous initial conditions when, suddenly, due 
to a rapid decrease of the temperature, the mass squared acquires a negative value. This model is an idealization of 
instabilities arising in, e.g., inflationary cosmology and during the chiral phase transition in expanding dense hadronic 
matter. Thus, we study the time evolution of the field described by the Hamiltonian 

H(t) = g(ftt(p)ft(p) + (m 2 (t) + p 2 ) &(p)$(pj) , (37) 

with 

2f+\ f to 2 > for t < , /00 N 

m (t) — < 2 ' (38) 

v ' \ -fi 2 < for t > . y ' 

At t < the state of the system is described by the equilibrium Wigner functional ( |l7| ) . To find the system evolution 
starting with t = 0, we split the Hamiltonian (B7j) for t > into stable and unstable modes as 



H(t > o) = £ ^(nt(p)n(p) - wl(p) &(p)Hp)] 

+ f ^(nt(p)n(p) + ul(p) &(p)Hp)) , (39) 



where to±(p) = yl ±(p 2 — fjp) is a real number. 

One sees that the modes are independent from each other, thus we discuss for a moment a single mode, stable or 
unstable. Since the Wigner functional equation of motion (Jl4|) with to 2 = — p, 2 coincides with the classical Liouville 
equation, it can be solved in a way which is well known from classical mechanics. Specifically, the single mode equation 
of motion is solved by 

W(%, U p ; t) = W($ p (-t), II p (-t); 0) , (40) 
where $ p (i) is the solution of the classical field equation 



G 



dt 2 



* p (f)±<4(p) $„(*) = 



(41) 



with n p (<) = &„(t) and the initial conditions $ p (0) = n p (0) = II p . 

Solving Eq. ([y]) one finds the single mode Wigner function (|io|), and then gets the whole functional as a product 
of the single mode functions. Thus, 



W[$,II;t] = Cexp 



where 



~0 J™ g a^(p) (ns(p)n„( P ) + £; 2 ( P )$5(p)$ (p) 



(42) 



$o(p) = ch(u;_(p)i) 



-(P) 



n(p) sh(w_ (p)t) 



n (p) = n(p) ch(w_(p)t) - cj-(p) $(p) sh(w_(p)i) 



for < p < /i and 



$o(p) = $(p) cos(co+(p)*) - 



w+(p) 



II(p) sin(w+ (p)t) 



IIo(p) = E(p) cos(w+(p)i) +w+(p) $(p) Bin(o; + (p)t) , 

for p > fi. It should be stressed that the solution is exact and fully quantum mechanical. 

If one is interested only in the dependence of the Wigner functional, II is integrated over and one gets 



VII 
~2T 



W[<f>,IL;t] = Cexp 



l Jq ZTT 



(43) 
(44) 

(45) 
(46) 



(47) 



where 



x(p,t) = Q(ji-p) 



+Q(p-H) 



w 2 _(p) ch 2 (w_(p)t) +E 2 (j>) sh 2 (uj-(p)t) 
E 2 {p)u>%(p) 



w\(p) cos 2 (u+{p)t) + E 2 {p) sin 2 (w+(p)t) 



Using Eq. 



one easily finds how the unstable modes grow 
(*t(p)$(p)) t = (*t(p)$(p)) fch 2 ( W _(p)t) 



(48) 



(49) 



with {ffi(p)$(p))o given by Eq. (pl|). The unstable mode growth is usually obtained directly from the classical 
equation of motion (El). In such a case however, one finds only the first term in the r.h.s of Eq. fl49|). The second 
term, which results from the proper incorporation of the initial condition and the interplay between $ and II in (0), 
is missing. Then it is argued that the zero-momentum mode grows the fastest. This statement, however, is not quite 
correct. Approximating the r.h.s of Eq. (|49|) as 

(&(p)Hp))t = (&( P )Hp))o (l + (M 2 + ™ 2 )t 2 + 0{uj 4 _(p)t 4 )) ■ (50) 

one sees that when the instability initially (for t <C u>Z (p)) all unstable modes (p 2 < p?) grow at the same rate, and 
the zero-momentum mode becomes dominant only at a later stage. This observation might be of physical significance 
since the experimental detection of the misaligned chiral condenstates [6], if they are produced in heavy-ion collisions, 
would be difficult if the phenomenon occurred only for pions with approximately zero momentum. Eq. (|5^) suggests 
that the pions with nonzero momentum have a chance to contribute to the condensate. 

The two-point correlation function is also of interest here, because it contains information about domain growth after 
the phase transition. The generating functional Z[j] defined in ( ps| ) depends explicitly on time in the present case. A 
straightforward calculation shows that the generating functional again has the form (pOj), where the time-dependent 
correlation function now is given by 
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1 ' J 2vr A^(p) X (p,t) V ' 

The late-time properties of Q(x,t) are determined by the exponentially growing modes with p 2 < [i 2 , for which 
x(p>t) ex P[ — 2o>— (p)*] ~^ 0. Expanding around p = and carrying out the three-dimensional Gaussian 

integral over p one finds: 



From this result one can read off that the size of correlated scalar field domains grows as (x 2 ) ~ At/ fi. 

VII. APPROXIMATE METHODS 



There are only a few situations when the exact equation of motion (11) can be solved analytically. Thus, there is 
an obvious need to develop approximate methods to study problems more complicated than those discussed above. 
We briefly present in this section three approaches which seem to be promising. 

A. Mean-field approximation 

The mean field approximation, which has been frequently used to study symmetry breaking phase transitions, is 
implemented by replacing in the initial Lagrangian the terms which are cubic, quartic, etc. in fields by the products 
of the fields and their expectation values. In the case of £/(<£>) = — |r <I> 4 the replacement is 

±fr(x)->j(&(x))&(x), (53) 

or equivalently, we replace the mass m 2 in the free field equation of motion by the effective mass m 2 = m 2 + -|($ 2 (x)). 
The combinatorial factor 4! changes into 4 in Eq. (B3) because there are 6 ways to select 2 fields out of 4. Since 
(^ 2 (x)) is determined by the functional W via Eq. (Bh, we get a self-consistent nonlinear equation for W. 

Let us briefly discuss the mean field approximation for the case of thermodynamical equilibrium. The Wigner 
functional is then given by Eq. (|l7]) derived for the free fields. However, instead of the mass there is the effective 
mass determined by the following gap equation 

ml = m 2 + ^ G{0) , (54) 

where £7(0) is defined by Eq. ( |3l] ) with to* substituting m. Before the equation ( |5^ ) is solved one has to regulate 
£7(0) subtracting from it a divergent zero temperature contribution. After this procedure we find the well known gap 
equation 



2 , A r dp ^Trnj 1 r ,_, 

m*=m +- — 7^7= =r o , o , (55) 



2 J-oo 27r exp(/Vp 2 + ml) - IP 

which is further analysed in |2] . It is worth mentioning that the standard result is reproduced only after the correlation 
function is renormalized. 

B. Semiclassical approximation 

This approximation might be particularly useful when the system of interest is, obviously, "semiclasical" and the 
"classical" solution, analytically or numerically, is known. Then, the quantum correction can be found in the following 
way. The equation of motion (|ll]) is written down in a symbolic notation as 

(k c + n 2 k q ) w = o, (56) 
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where K c and K q represent the "classical" and "quantum" part of the operator acting on the Wigner functional, 
which is expanded as W = W c + h 2 W q . Since K C W C = 0, we find W q as a solution of the equation 

K c W q = -K q W c + 0(h 2 ) . (57) 
For a single particle Wigner function one finds 

W q (q,p;t) = J dt'dq'dp 1 G(q,q',p,p';t,t') K q W c (q' ,p';t') (58) 

where the Green function G(q, q' ,p,p']t,t') solves the equation 

K c G(q, q',p,p'; t, f ) = — + — — - — — = 6(t - if) 5{q - q') 5{p - p') . (59) 

Since solutions of the classical equations of motion are assumed to be known, the Green function can be explicitely 
written as 

G(q,q',p,p';t,t') = Q(t - if) S(q - q(t) - q' + q{if)) 5[p - p{t) - p' + p{t')) , (60) 

with (q(t),p(t)) being the classical solution. The generalization of Eqs. (|58|) and ( |60[ ) to the case of a quantum field 
is trivial, only the notation is more complicated. 



C. Numerical Simulation of the Initial Value Problem 



The main goal of this study has been to find a formalism to analyze temporal evolution of the quantum fields 
beyond the classical limit. Here we argue that equation of motion ([ll]) is indeed numerically tractable. 

The first step is to define the field $(x) on a lattice. Then, the equation of motion ( pi] ) resembles a many particle 
transport equation with every site of the lattice representing a single particle that interacts with an external polynomial 
potential (due to the mass and interaction terms of the Lagrangian) and with its nearest neighbors (via the gradient 
term). We can then adopt, e.g., the method developed by John and Remler to solve the quantum transport 
or quantum Liouville equation of a many particle system. The method is sketched below for a single particle; a 
generalization to many particles is straightforward. 

A single particle Wigner function is represented by a finite sum over discrete points in phase space. Since the Wigner 
function is not positive definite, the points may have positive or negative signature. Then, each point representing 
a test particle is evolved in such a way that every step of the "classical" evolution is followed by "quantum" step, 
which converts a single momentum point into several ones. Specifically, the delta function in momentum space is 
represented by a gaussian distribution centered around the classical position. The higher momentum gradients, which 
are responsible for the quantum evolution, act on this distribution. The resulting function, which has positive and 
negative components, is again represented by discrete points, i.e. by new test particles with positive or negative 
signature. To avoid unlimited proliferation of the test particles, those with closely neighboring trajectories can be 
allowed to merge. It was demonstrated in [ pcfl that the method successfully works provided the initial Wigner function 
is sufficiently smooth in phase space. We intend to apply this method in future work to study the quantum field 
evolution. 



VIII. SUMMARY AND OUTLOOK 



We have introduced the formalism which describes the quantum field dynamics in terms of the field amplitudes $ 
and their conjugate momenta EE. The Wigner functional playing the central role in our approach provides a density in 
a phase space spanned by $ and IT. We have derived the equation of motion of the functional, which is of the form of 
the quantum transport or Liouville equation. It has been shown in the last section how the knowledge of the classical 
solution of this equation can be used to find quantum corrections. The free fields in thermodynamical equilibrium 
have been discussed in detail, and we found an explicit analytical solution of the field analog of the "upside-down" 
harmonic oscillator. The solution incorporates in a very natural way the "thermal" initial conditions, its physical 
meaning is intuitive and transparent. In the last section we briefly discussed the approximate methods which will be 
needed to solve more complicated problems. 

We intend to apply our approach to study the role of quantum fluctuations in temporal evolution of the symmetry 
breaking phase transitions beyond the gaussian approximation [Q. As mentioned in the introduction our formalism 
also appears useful to investigate the quantum content of solutions of the classical field equations of motion. 
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APPENDIX A: 



In this Appendix we prove the equations (||) and (^), which show that the Wigner functional can be treated as a 
density in the phase space spanned by <& and II. The formula (0) is obvious if one observes that 



exp 



i / dx U(x)(/)' (x) 



= 8{4>'{x)) 



(Al) 



When the operator O depends only of $ but not ft, one easily gets Eq. (g) since the states |$) are by definition 
eigenstates of $ and consequently 0($)|$) = 0(<fr)\<&). It is assumed that C($) can be expanded in the power series 
of l>. 

Let us now consider O which depends only on II. Then, we introduce the complete set of the momentum eigenstates 
and the r.h.s of Eq. (iq) is 



i r T)TT T)T\ T)T\ f 

J_ J V^Vcp'— 2?r 1 2 J 0(IL 2 )exp[-ij dxU{x)<t>i{x)] 



x ($ + i^'iniXHii p(t) |n 2 )(n 2 |$ - \<j>') 



(A2) 



Further we use the momentum eigenfunction 

($|n) = exp 

and observe that 

0(U 2 (x)) exp 



i / dx H(x)$(x) 















exp 


ijdx( 


n x -n a 


*4 


ni + n 2 


*)1 



oil 4 



1 



iS(j)'(x) 2i5^{x) i 















) exp 


ijdx( 


ni-n 2 


*4 


Hi + n 2 


*)] 



(A3) 



(A4) 



Performing repeatedly the partial integration we effectively convert Oyj g^f^ — ji sJ( x ) J Ofa( x )) an d finally 

prove Eq. (§) for 0(lL(x)). 

When the operator O depends simultaneously on $ and II, the situation gets somewhat complicated because the 
operators $ and II do not commute with each other. The direct computation, in particular, shows that 



(TL(xMy)) = (Il(x)$(y)) - - 5(x - y) , 



(A5) 



while 



Therefore, 



where {A, B} = h [AB + BA 



($(y)Il(x))=(Il(x)$(y)) + -5(x-y) 



({%),<%)}) Hn(x)<%)) , 



(A6) 
(A7) 



Generalizing the result (A.7), one proves the equality (C(<f>,II)) = (C(<f>, II)) assuming that the pairs of noncom- 
muting operators are symmetrized. 
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APPENDIX B: 



We derive here the equation of motion (3.3). For this purpose the Hamiltonian ( [To| ) is split as 

H = Hn + Hy + H m + Hi , 



(Bl) 



where Hjj, H^, H m , and Hj correspond to the first, second, third and fourth term, respectively, on r.h.s of Eq. jl^ ) 
Now one has to calculate four expressions 



Gi = / V(j>' exp 



-- / dxU(x)(f)'(x) 



with the index a — II, V, m and /. 

The evaluation of G m is straightforward since |<I>) is, by definition, the eigenstate of $. Thus, 



G m = — I dx I V(j)' exp 



•- / dx n(x)^'(x) 



x (($(*) + l -ct>\x)f - ($(x) - y\x)) 2 ) ($ + y\ p i$ - y) 



Since, 



one finds 



4>'(x) exp 



dx U(x)(j)'(x) 



ih- 



SU(x) 



exp 



dx U(x)(p'(x) 



G m = ihm 2 [ dx $(x) —^T W[$,Il;i\ 
J 6Tl(x) 



Keeping in mind that Ci is a polynomial in $, one also easily computes Gi as 



G, 



dx LC/($(x) 



ih s 



-£i($(x)- 



ih 6 



2 SU(x)J ~ J 1*" 2 SU(x) 
To find Gn we introduce the complete set of the momentum eigenstates |II) and then 

dil vn 2 



w[*,n-,t] 



Cn =- IdxJ — — ex p 

1 .,.„ \ /„ , 1 



l - J dx n(x)<j>' (x) 



{Uj(x)-U 2 2 (x)) 



x (* + -0 / |niHni|p|n 3 )(n a |$--0 / 



Next we make use of the explicit form of the momentum eigenf unctions (A3) and observe that 



27T 



i 2 6$(x) 8<f>'(x) 



exp 



n 1 lLr 



(n a -n 2 )$ + -(ni + n 2 y 



= (nf - n£) exp 



(n 1 -n 2 )$ + ~(n 1 + n 2 y 



(B2) 



(B3) 



(B4) 



(B5) 



(B6) 



(B7) 



(B8) 



Since the derivative over <fi' appears under the integral over <fi' , we perform the partial integration and as a result we 
obtain 



G n = -ih [ dx U(x) -^-r W[$, II; t] . 
J d<P(x) 

We best find Gy defining the field <i>(x) on a lattice. The Hamiltonian is then 



(B9) 



a x -\ 

9 2^ 



(BIO) 
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where a is the lattice spacing and i numerates the lattice sites. With such a Hamiltonian the computation of Gv is 
very similar to that one of G m and one finds 



g v = ma £ - * 4) (JL- - A) w(m , m . t) 

In the continuum limit the above expression reads 

v$(x) v^^y n ; i] . 

After the partial integration we finally get 

g v = -in f dx vH( x ) w[&, n ; t] . 

J 5Tl{x) 
Collecting Gn , Gv , G rn and Gj we find the desired equation of motion (fill) . 



(Bll) 



(B12) 



(B13) 



APPENDIX C: 

We derive here the equilibrium Wigner functional of the fields in the momentum space. We introduce the complete 
set of the energy eigenstates of the two-dimensional (isotropic) oscillator \ni,ri2) and rewrite Eq. (0) as 



EI] = / V(t>' exp 



x (& + ~<t>'\n 1 ,n 2 ){n2,n 1 \§- ^<t>') 



(CI) 



with E ni _, l2 (p) = \/p 2 + m 2 (m + ri2 + 1). The eigenfunctions for the two dimensional (isotropic) oscillator are |21 

~,2 \ 1/2 



($ fl ,$ z ni,n 2 ) 



7r 2"i n x ! 2™2 n 2 



x H ni (a$ R ) H n2 {a$i) exp 



where $^,7 is the real or imaginary part of the field, respectively. H n denotes the Hermite polynomial, while 



9 2 + m 2 \ V4 



Substituting ( |C2| ) into ( |Cl|) and using the identity |19| 

°° a n 1 
V — H n (x)H n (y) = exp 
^ VI - 4a 2 



Aaxy - 4t 2 (x 2 + y 2 ) 
1 - 4a 2 



(C2) 



(C3) 



(C4) 



which holds for a < 1/2, one finds after elementary integration and proper normalization the final formula of the 
equilibrium Wigner functional @. The same result can be obtained by analytic continuation and generalization to 
infinitely many degrees of freedom of the formula given in the Appendix of ref. [l^ ] for the generating function G(s) 
for the Wigner functions associated with the energy eigenstates of the harmonic oscillator. 



APPENDIX D: 

The purpose of this Appendix is to derive the generating functional (|3^). We apply the method described in the 
literature [22| modifying it slightly. After integration over II and partial integrations with respect to x' and then x, 
we rewrite Eq. (p8|) as 
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Z[j] = C / X>$ exp 



-0 I dx dx'[-A(x - x')<f>(x)(m 2 - V 2 )${x') - 6(x - x')$(x)j(x') 



(Dl) 



It is important for these manipulations that A(x) = A(— x). Now we change the variable into <I>(x) + $o(x) 

demanding that the held <&o(x) satishes the equation 

(D2) 



(D3) 



dx'A(x - x 1 ) ( - V 2 + m 2 ) $ Q (x') = j(x) . 
After integration over $ and the partial integrations with respect to x and x' ', we get the result 



Z [j] = W exp 



| / </-<■ <!>< i (■'■)./(.'■) 



Substituting into Eq. (D3) the solution of Eq. (D2) in the form 

$o(z) = / dx' G(x - x')j(x') , 
with Q being the Green function given by Eq. (|3l|), one hnally finds the generating functional ( (30| 



(D4) 
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